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Abstract. We consider a long-range Ising antiferromagnet (LRIAF) put in a transverse field. 
Applying quantum Monte Carlo method, we study the variation of order parameter (spin 
correlation in Trotter time direction), susceptibility and average energy of the system for 
various values of the transverse field at different temperatures. The antiferromagnetic order 
is seen to get immediately broken as soon as the thermal or quantum fluctuations are added. 
We also discuss the phase diagram for the Sherrington-Kirkpatrick (SK) model with the same 
LRIAF bias, also in presence of a transverse field. We find that while the antiferromagnetic 
order is immediately broken as one adds an infinitesimal transverse field or thermal fluctuation 
to the system, an infinitesimal SK spin glass disorder is enough to induce a stable glass 
order in the antiferromagnet. This glass order eventually gets destroyed as the thermal or 
quantum fluctuations increased beyond their threshold values and the transition to para phase 
occurs. Indications of this novel phase transition are discussed. Because of the presence of full 
frustration, this surrogate property of the LRIAF for incubation of stable spin glass phase in 
it (induced by addition of a small disorder) should enable eventually the study of classical and 
quantum spin glass phases by using some perturbation theory with respect to the disorder. 



1. Introduction 

Quantum phases in frustrated systems are being intensively investigated these days; in particular 
in the context of quantum spin glass and quantum ANNNI models pp. Here we study a fully- 
frustrated quantum antiferromagnetic model. Specifically, the long-range antiferromagnetic Ising 
model put under transverse field. The finite temperature properties of sub-lattice decomposed 



version of this model was already considered earlier [3j 0] . The quantum phase transition and 
entanglement properties of the full long-range model at zero temperature was studied by Vidal et 
al [5]. Here we present some results obtained by applying quantum Monte Carlo technique [6] to 
the same full long-range model at finite temperature. We observe indications of a quantum phase 
transition in the model, where the antiferromagnetically ordered phase gets destabilised by both 
infinitesimal thermal (classical) as well quantum fluctuations (due to tunneling or transverse 
field) and the system becomes disordered or goes over to the para phase. 

The ordered phase of the long range Ising Antiferromagnet (LRIAF) seems to be extremely 
volatile and loses the order (freezing of spin orientations) at any finite fluctuation level; classical 
or quantum. However the LRIAF model has the required frustration of the Sherrington- 
Kirkpatrick (SK) model, which could support the spin glass order, but for any disorder. To 
check if this 'liquid'-like antiferromagnetic phase of LRIAF can get 'crystallized' into spin-glass 
phase if a little disorder is added, we study next this LRIAF Hamiltonian with a tunable coupling 
with the SK spin glass Hamiltonian and study this entire system's phase transitions induced by 
both thermal and tunneling field. Indeed, the stable SK spin glass phase is observed for thermal 
or quantum fluctuations below finite threshold values. 

We employ the analytic (mean field) solution of the transverse Ising model with long-range 
interactions can be in presence of the special kind of quenched disorder appropriate for the SK 
spin glasses. With this, we study the phase diagram for the SK model with LRIAF bias in a 
transverse field. We again find that the antiferromagnetic order is immediately broken when 
one adds an infinitesimal transverse field or thermal fluctuation to the system. However, an 
infinitesimal SK-type disorder is enough to make the system 'crystallized' into a glass phase. 

This paper is organized in the following manner. In Section 2, we introduce the pure quantum 
LRIAF model and then discuss the (finite temperature) quantum Monte Carlo results. In Section 
3, we consider the SK model with antiferromagnetic bias in a transverse field, and discuss the 
(analytic) mean field phase diagram. In Section 4, we present some discussions on our results. 

2. The pure LRIAF model 

The Hamiltonian of the infinite-range quantum Ising antiferromagnet (without any spin glass 
disorder) is 

T N N N 

h ee n<n+n<n = ± £ ^-/^f -r^>f, (i) 

ij'(>i)=l i=l i=l 

where Jo denotes the long-range antiferromagnetic (Jo > 0) exchange constant; for convenience, 
we fix the value Jo = 1 in this section. Here a x and a z denote the x and z component of the N 
Pauli spins 

<=(J °f = (? J); < = i,2,....,iv. 

h and T denote respectively the longitudinal and transverse fields. We have denoted the co- 
operative term of H (including the external longitudinal field term) by H^ c > and the transverse 
field part as . As such the model has a fully frustrated (infinite-range or infinite dimensional) 
co-operative term. At zero temperature and at zero longitudinal and transverse fields, the 
H( C ' would prefer the spins to orient in ±z directions only with zero net magnetization in 
the z-direction. This antiferromagnetically ordered state is completely frustrated and highly 
degenerate. Switching on the transverse field T would immediately induce all the spins to orient 
in the x-direction (losing the degeneracy), corresponding to a maximum of the kinetic energy 
term and this discontinuous transition to the para phase occurs at T = 0. However, at any finite 




Figure 1. Variation of the order parameter r (correlation in the Trotter direction) with transverse field V for 
T = 0.10,0.20 and 0.30 (h = 0) for two different system sizes (TV = 100 and 200). r = for large V. The inset 
shows the plot of r against the scaled variable F/T. 
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Figure 2. Variation of the susceptibility \ with transverse field Y for T = 0.10,0.20 and 0.30 (h < 0.1) for 
two different system sizes (TV = 100 and 200). The corresponding susceptibility \ci for various temperatures for 
TV = 100 and 200 for the classical system are shown in the inset. \ converges to the classical values \ci for large 

r. 



temperature the entropy term coming from the extreme degeneracy of the antiferromagnetically 
ordered state and the close-by excited states does not seem to induce a stability of this phase. 

2.1. Monte Carlo simulation 

2.1.1. Suzuki- Trotter mapping and simulation This Hamiltonian ([T]) can be mapped to a 
(oo + l)-dimensional classical Hamiltonian [6] using the Suzuki- Trotter formula. The effective 
Hamiltonian can be written as 

1 N p , TV p j N P 

n = Yj Yl a ^i,k - p Yl ai ' k - ~p Yl Yj a ^i,k+i , 

ij(>i)=l k=l i=l fc=l i=l k=l 
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Figure 3. Variation of average energy E with transverse field r for T = 0.10, 0.20 and 0.30 (h = 0) for two 
different values of N(— 100,200). The corresponding average energy E c i for various temperatures for TV = 100 
and 200 for the are shown in the inset. E converges to the classical values Ed for large F . 



where 

J p = -(PT/2) ln(tanh(r/Pr)). (2) 

Here P is the number of Trotter replicas and k denotes the fc-th row in the Trotter direction. J p 
denotes the nearest-neighbor interaction strength along the Trotter direction. We have studied 
the system for N = 100. Because of the diverging growth of interaction J p for very low values 
of r and also for high values of P, and the consequent non-ergodicity (the system relaxes to 
different states for identical thermal and quantum parameters, due to frustrations, starting from 
different initial configurations), we have kept the value of P at a fixed value of 5. This choice 
of P value helped satisfying the ergodicity of the system up to very low values of the transverse 
field at the different temperatures considered T = 0.10 and 0.20. Starting from a random initial 
configurations (including all up or 50-50 up-down configurations) we follow the time variations 
of different quantities until they relax and study the various quantities after they relax. 

2.1.2. Results We studied results for three different temperatures T = 0.10,0.20 and 0.30 and 
all the results are for N = 100 and 200 and P = 5. We estimated the following quantities after 
relaxation : 

(i) Correlation along Trotter direction (r) : We studied the variation of the order parameter 

. N P 

r = J^^Z^Z^i^k+l) , (3) 
t=l fc=l 

which is the first neighbor correlation along Trotter direction. Here, (...) indicate the average 
over initial spin configurations. This quantity r shows a smooth vanishing behavior. We consider 
this correlation q as the order parameter for the transition at r c . A larger transverse field is 
needed for the vanishing of the order parameter for larger temperature. The observed values 
(see Figure H]) of T c are ~ 1.6,2.2 and 3.0 for T = 0.1,0.2 and 0.3 respectively. As shown in 
the inset, an unique data collapse occurs when r is plotted against T/T and one seems to get 
the complete disorder immediately as the scaling dos not involve any finite value T c . This is 
consistent with the observations in the next section. 



(ii) Susceptibility (\) ■ The longitudinal susceptibility x = (V^-P^Ei k( a hk)]/9h, where 
h (— > 0) is the applied longitudinal field, has also been measured. We went up to h = 0.1 and 
estimated the x values. As we increase the value of the transverse field T from a suitably chosen 
low value, x initially starts with a value almost equal to unity and then gradually saturates at 
lower values (corresponding to the classical system where J p = in Eq.(2)) as T is increased. 
Also at r = 0, the classical values are indicated in Figure [21 This saturation value of x decreases 
with temperature. Again the field at which the susceptibility saturates are the same as for the 
vanishing of the order parameter for each temperature. 

(iii) Average energy (E) : We have measured the value of the co-operative energy for each 
Trotter index and then take its average E i.e. E = (H^) of Eq. (1) with J = 1. It initially 
begins with —1.0 and after a sharp rise the average energy saturates, at large values of T, to 
values corresponding to the classical equilibrium energy {E c i for J p = in Eq.(2)) at those 
temperatures. Again it takes larger values of T at higher temperatures to achieve the classical 
equilibrium energy. At T = 0, the corresponding classical values of E are plotted in Figure 
O The variations of all these quantities indicate that the 'quantum order' disappears and the 
quantities reduce to their classical values (corresponding to J p = in for large values of the 
transverse field T. 



2.2. A mean field analysis 

First, Let us consider the case of pure LRIAF model and rewrite our Hamiltonian H in Eq.(l) 
for J = as 



/ N \ z N N N 

\?:=i / i=i i=i i=i 



(4) 



If we now denote the total spin by Btot i-e. fftot = j? Y2i=i &i (where N\a\ = 0, 1, 2, N), then 
the Hamiltonian H can be expressed as 

-jy = -^tot) 2 ~ hcrl ot - Taf ot - — . (5) 

Let us assume the average total spin (a) to be oriented at an angle 9 with the z-direction : 
( a tot) = mcosO and (af ot ) = msm.9. Hence the average total energy E to t = (H) can be written 

as El 1 

— — = —m 2 cos 2 6 — hmcos 9 — Tm sin 9 . (6) 

N 2 N K ' 

At the zero temperature and at V = 0, for h = 0, the energy Etot is minimised when 9 = 
and m = (complete antiferromagnetic order in z-direction). As soon as T ^ (h = 0) the 
minimisation of Etot requires 9 = n/2 and m = 1 (the maximum possible value); driving the 
system to paramagnetic phase. This discontinuous transition at T = was also seen in [5]. 
As observed in our Monte Carlo study in the previous section, T C (T) — > as T — > 0. This 
is consistent with this exact result r c = at T = 0. For T = (and h = 0), therefore, the 
transition from antiferromagnetic (0 = = m) to para {9 = 7r/2,m = 1) phase, driven by the 
transverse field T, occurs at T = itself. 

One can also estimate the susceptibility x at T = = T. Here Etot IN = ^m 2 cos 2 # — 
hmcos 9 — -k and the minimisation of this energy gives m cos 9 = h giving the (longitudinal) 
susceptibility x = mcos9/h = 1. This is consistent with the observed behaviour of x shown in 
Figure^] where the extrapolated value of x a t T = increases with decreasing T and approaches 
X = 1 as T -> 0. 

At finite temperatures T / 0, for h = 0, we have to consider also the entropy term and 
minimise the free energy T = Etot — TS rather than Etot where S denotes the entropy of the 



state. This entropy term will also take part in fixing the value of 9 and m at which the free 
energy J- is minimised. As soon as the temperature T becomes non-zero, the extensive entropy 
of the system for antiferromagnetically ordered state with m ~ (around and close-by excited 
states with 9 = 0) helps stabilisation near 9 = and m = rather than near the para phase 
with 9 = 7r/2 and m = 1, where the entropy drops to zero. While the transverse field tends to 
align the spins along x direction (inducing 9 = ir/2 and m = 1), the entropy factor prohibits 
that and the system adjusts 9 and m values accordingly and they do not take the disordered 
or para state values [9 = n/2 and m = 1) for any non-zero value of V (like at T = 0). For 
very large values of T, of course, the free energy T is practically dominated by the transverse 
field term in H and again 9 = n/2 and m = 1, beyond T = T C (T) > for T > 0. However, 
this continuous transition-like behaviour may be argued [7J to correspond to a crossover type 
property of the model at finite temperatures (suggesting that the observed finite values of T C (T) 
are only effective numerical values). In fact, for h = one adds the entropy term — TlnD s to 
E to t in Eq.(7) to get T and one can then get [7J, after minimising the T with respect to m and 
9, m = tanh(r/2T), which indicates an analytic variation of m and no phase transition at any 
finite temperature for J = (antiferromagnetic phase occurs only at T = T = as shown in 
Figure [U 



3. LRIAF with SK disorder: 'Liquid' phase of the SK spin glasses gets frozen 

In the last part of this paper, we discuss the phase diagram for the Sherrington-Kirkpatrick model 
with antiferromagnetic bias in a transverse field. We find that the antiferromagnetic order is 
immediately broken when one adds an infinitesimal transverse field or thermal fluctuation to 
the system, whereas an infinitesimal SK-type disorder is enough to get the system 'crystalized' 
into the glass phase. 

The model we discuss here is given by the following Hamiltonian 

H = jf E (-W^K^-r^of (7) 

ij(j>i) i 

where Jo is a parameter whcih controlls the strength of the antiferromagnetic bias and J 
is an amplitude of the disorder Ty in each pair interaction. The T controls the quantum- 
mechanical fluctuation. When we assume that the disorder Tjj obeys a Gaussian with mean 
zero and variance unity, the new variable Jjj = — Jo + Jr^- follows the following distribution. 
P{Jij) = exp[-(J i + J ) 2 /2J 2 ]/V2ttJ. Therefore, we obtain the 'pure' antiferromagnetic Ising 
model with infinite range interactions when we consider the limit J — > keeping Jo > 0. On 
ther other hand, for Jo < with T = is identical to the classical SK model. In this paper, we 
investigate the condition for which the antiferromagnetic order survives. 

For the Hamiltonian d7J), we immediately obtain the saddle point equations under the static 
and replica symmetric approximations as follows (see e.g. [21 E] ) - 



m 



I" Dz ( J VQ Z + K ° m ) tanh pJtf^z + KomF + T* (8) 
^{J^qz + Komy + T 2 

2 



/oo 
-oo 



Dz{ - ( J v^ + Er o m ) V, tanh 2 /3 A /(«/^ + ^om) 2 + r2, (9) 
( J^qz + K m) 2 + T 2 



where Kq = sgn(Jo)|Jo| and in = N 1 Yli a i ^ s a magnetization and q = N 1 Ylii l \ a t) 2 ls a 
spin glass order parameter. We defined Dz = dze~ z / 2 /\/2~k. The bracket (• • • ) denotes an 



expectation over the density matrix: p = e~@ H /tr e~^ H . When Jo is negative, ([8]) has the only 
solution m = 0. The general phase boundaries (see Figure 4) between the ferro (F), spin glass 
(SG), antiferro (AF) and para (P) phases can be obtained by solving the the above two equations 
in the limit m — > o; q ^ for the F-SG boundary, q — > for SG-P boundary and m — > fot the 
F-P boundary. In these limits, the P-SG boundary equations become (see e.g. [8]) 

r = Jtanh 

3.1. Classical system 

In the classical limit, the equations of state are simplified as m = Dztanh f3(J^/qz — Jom), 
q = D z tanh 2 (i{J^/qz - J m). For J > 0, we find that m = is only physical solution 
for all temperature regimes. This fact means that there are three possible phases, namely, 
the antiferromagnetic phase, the paramagnetic phase and the spin glass phase. In these three 
phases, the magnetization m is zero. To determine the critical point T$g at which the spin 
glass transition takes place, we expand the equation with respect to q for q ~ and m = 0. 
Then, we have T$g = J and the critical point is independent of the antiferromagnetic bias J> 





Figure 4. The phase diagram of classical SK model [5] extended for antiferromagnetic bias. For Jo < 0, 
there exist spin glass phase below T/J = 1 and the critical temperature is independent of the strength of the 
antiferromagnetic bias Jo. For finite temperature T > 0, the anti-ferromagnetic order dissapears and the system 
changes to the paramagnetic phase. When we add an infinitesimal disorder J > 0, the antiferromagnetic order is 
broken down and the system suddenly gets 'crystalized' into a spin glass (SG) phase. 

This result means that the antiferromagnetic order can apprear if and only if we set Jo > and 
T/J = 0, Jo/ J = 0. On the other hand, for < Jo < oo at low temperature regime T < T$g, 
the spin glass phase appears. We plot the phase diagram in Figure [H We also conclude that 
the system described by the Hamiltonian ([7]) with T = is immediately 'crysterlized' when we 
add any infinitesimal disorder J > 0. 

From the view point of the degeneracy of the spin configurations, we easily estimate the 
number of solution for the antiferromagnetic phase as 2 N ^ 2 = e °- 346Ar , which is larger than the 
number of the SK model e°- 199N [gj. However, for the infinite range antiferromagnetic model, 
the energy barrier between arbitrary configurations which gives the same lowest energy states 
is of order 1 and there is no ergodicity breaking. 



3.2. Quantum system 

We next consider the case of presence of transverse field r / 0. In this case, we also find that 
the saddle point equation ([8]) has a solution m = and the phase boundary between the spin 
glass and paramagnetic phases is given by setting m = and q ~ and we get (fTOj) . Obviously, 
the boundary at T = gives T$g = J ■ On the other hand, when we consider the case of V ~ 0, 
we have T$g = J- These fact means that there is no antiferromagnetic nor the spin glass phase 
when we consider the pure case J = because the critical point leads to T$g = Tsg = 0. 
Therefore, we conclude that the antiferromagnetic phase can exisit if and only if T = V = 0. 




Figure 5. Phase diagram for the quantum system. The antiferromagnetic order exisits if and only if we set 
T = r = 0. As the J decreases, the spin glass phase gradually shrinks to zero and eventually ends up at an 
antiferromagnetic phase at its vertex (for V = — T — J) as discussed in Section 2. 



4. Discussion 

We considered here a long-range Ising antiferromagnet at a finite temperature and put in a 
transverse field. The antiferromagnetic order is seen to get immediately broken as soon as the 
thermal or quantum fluctuations are added. However, when we add the Sherrington-Kirkpatrick 
hamiltonian as perturbation we find that an infinitesimal SK spin glass disorder is enough to 
induce a stable glass order in this LRIAF antferromagnet. This glass order eventually gets 
destroyed as the thermal or quantum fluctuations increased beyond their threshold values 
and the transition to para phase occurs. As shown in the phase diagram in Figure 5, the 
antiferromagnetic phase of the LRIAF (ocurring only at J = = T = T), can get 'crystalised' 
into spin-glas phase if a little SK-type disoder is added (J / 0); the only missing element in the 
LRIAF (which is fully frustrated, but lacks disorder) to induce stable order (freezing of random 
soin orientations) in it. As mentioned already, the degeneracy factor e 0,346Ar of the ground state 
of the LRIAF is much larger than that e 0,199Ar for the SK model. Hence, (because of the presence 
of full frustration) the LRIAF posses a surrogate incubation property of stable spin glass phase 
in it when induced by addition of a small disorder. This observation should enable eventually 
the study of classical and quantum spin glass phases by using some perturbation theory with 
respect to the disorder. 
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